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SUMMARY 


A procedure is outlined for the numerical solution of the complete elas- 
tohydrodynamic lubrication of rectangular contacts incorporating a non- 
Newtonian fluid model. The approach uses a Newtonian model as long as the 
shear stress is less than a limiting shear stress. If the shear stress exceeds 
the limiting value, the shear stress is set equal to the limiting value. The 
numerical solution requires the coupled solution of the pressure, film shape, 
and fluid rheology equations from the inlet to the outlet. Isothermal and 
no-side-leakage assumptions were imposed in the analysis. 

The influence of dimensionless speed U, load W, materials G, and 
sliding velocity U* and limiting-shear-strength proportionality constant y 
on dimensionless minimum film thickness H m ^ was investigated. Fourteen 
cases were used in obtaining the minimum-film-thickness equation for an elas- 
tohydrodynamical ly lubricated rectangular contact incorporating a non-Newtonian 
fluid model 
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where 


\ln # N " 3,07 U 0 * 71 G 0,5 V°* 11 

Computer plots are also presented that Indicate In detail pressure dlstrlbu 
tlon, film shape, shear stress at the surfaces, and flow throughout the 
conjunction. 


INTRODUCTION 

To obtain a better understanding of the failure mechanism In machine 
elements, the next generation of elastohydrodynamlc lubrication analysis 
should Incorporate such effects as 

(a) Non-Newtonian fluid 

(b) Surface roughness 

(c) Temperature 

Although a Newtonian solution for elliptical contacts had been obtained by 
Hamrock and Dowson (1981), it was felt that a rectangular or line-contact 
analysis should be used to incorporate these effects because of their added 
complexity. A rectangular or line-contact analysis was performed by Hamrock 
and Jacobson (1982). The analysis required the simultaneous solution of the 
elasticity and Reynolds equations. The equations were coupled from the inlet 
to the outlet assuming Isothermal conditions and no side leakage. The re- 
sults from Hamrock and Jacobson (1982) are the foundation for the more com- 
plicated analysis Incorporating non-Newtonian fluid, surface roughness, and 
temperature. When considering these effects, the Initial pressure profiles 
are those obtained from the authors' earlier work. 

The present paper considers incorporating a non-Newtonian fluid model 
Into the theory of elastohydrodynamlc lubrication of rectangular contacts. 

The geometry of the problem is a roller, rolling and sliding against a plate. 
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where tht roller length it Urge as compared with tht radius. Tht lubricant 
In in elastohydrodynamlc conjunction experiences rapid and vtry large pres- 
turt variations, • rapid transit tint, possible largo temperature changes, and 
particularly In sliding contacts, high shoar rates. Tht groat severity of 
thoso conditions has cal ltd Into question tho normal assumption of New- tonlan 
behavior. The approach to be used In this paper Is to redefine the pressure 
and mass flow rate equations depending on how the values of shear stress at 
the surfaces compare with the limiting shear stress. The limiting shear 
expression used Is a semlemplrlcal linear function of pressure. 

Gaclm and Winer (1981) used a non-Newtonian fluid rheological model In 
their elastohydrodynamlc lubrication studies. Some limitations of this work 
are listed below: 

(1) The analysis assumes a non-Newtonian fluid model for the entire 
conjuctlon Including the Inlet region. The present paper uses a non- 
Newtonian fluid model only when the shear stress at the surfaces exceeds the 
limiting shear stress. 

(2) The Gaclm and Winer analysis relies on using a Grubln type of 
solution rather than employing a complete solution of the Reynolds, rheology, 
and film shape equations as Is used In the present paper. 

(3) The Gaclm and Winer paper assumes that the limiting shear stress Is 
zero when the fluid pressure Is zero. The present paper assumes that the 
limiting shear stress Is equal to an Initial shear strength when the pressure 
Is zero. 

figure 1 shows the effect of shear stress on shear strain rate for the 
present model and that of a Newtonian fluid. From this figure It Is observed 
that, In the present model, If the Newtonian shear stress exceeds the limit- 
ing shear stress, the shear stress Is set equal to the limiting shear stress. 
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The fluid model It Newtonian except when the sheer stress retches the sheer 
strength velue. At thet point, slippage occurs end the sheer stress Is equel 
to the sheer strength. 

Besides the dimensionless loed, speed, end meterlels perimeters thet were 
found to heve en Influence on film thickness In the authors' eerller peper 
(Hamrock end Jtcobson, 1982), when non-Newtonian effects ere considered, two 
eddltlonel peremeters were found to Influence the minimum film thickness, 
nemely: 

(1) Sliding velocity 

(2) Limiting-shear-strength proportionality constant 

Fourteen cases were used In obtaining a fully flooded film thickness 
equation when considering non-Newtonian effects of the liquid. Besides the 
film thickness calculations that were made, calculations of the force com- 
ponents, shear forces, coefficient of friction, and center of pressure were 
also performed. Computer plots are presented that Indicate pressure distri- 
bution, lubricant film shape, flow, and shear stresses within the conjunction. 

SYMBOLS 

B 1/n 

b semiwidth of Hertzian contact, rVbw/ 7, m 

¥ b/n, m 

E modulus of elasticity, N/m £ 

E' effective elastic modulus, 2/[(l - v 2 )/E t * (1 - v|)/ E b ], N/m 2 

7 dimensionless shear force 

f shear force per unit length, N/m 

G dimensionless materials parameter, «E' 

H dimensionless film thickness, h/R 

dimensionless minimum film thickness, h m ^ n /R 
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TL< n dimensionless minimum film thickness obtained from 
least-squares fit of data 

'tyirin h dimensionless minimum film thickness obtained from 

* least-squares fit of data while assuming a Newtonian fluid 
model 

h film thickness, m 

h^ minimum film thickness, n: 

m dimensionless distance from inlet to the center of Hertzian 

contact 

n number of nodes In semi axis of contact 

P dimensionless pressure, p/E* 

2 

p pressure, N/nr 

Q dimensionless mass flow per unit length, o/pqU $ R 

q mass flow per unit length, kg/(s«m) 

R effective radius In x direction, m 

r curvature radius, m 

S geometrical separation, m 

U dlmenslonler i speed parameter, ngU s /E'R 

U dimensionless velocity, u/uj 

U* dimensionless sliding velocity, uJu e 

a s ORIGINAL PAGE IS 

u velocity In direction of motion, m/s OF POOR QUALITY 

Uj velocity difference, (u # - u b )/2, m/s 

u s velocity sum, (u # + u b )/2, m/s 

Yi (Hnln - ^1n)100/4ui n 

V3 (p - p) 100/p 

W dimensionless load parameter, w z /E'R 

load per unit length, N/m 

X dimensionless coordinate, x/b 
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X 

Z 

z 

a 

Y 
6 
c 
n 
n 

n 0 

u 

V 
p 
"p 

p 0 

t 

T 

T 

t l 

t l 

T o 

Subscripts: 

a 

b 

x 

z 


coordinate in direction of motion, m 

* OF POOR 

dimensionless coordinate In direction of film, z/h 
coordinate In direction of film thickness, m 
pressure-viscosity coefficient of lubricant, m 2 /N 
limiting-shear-strength proportionality constant 
elastic deformation, m 
coefficient of determination 
absolute viscosity at gage pressure, N s/m? 
dimensionless viscosity, ti/hq 


PAGE 19 
QUALITY 


viscosity at atmospheric pressure, N s/m 2 
coefficient of friction 
Poisson's ratio 
lubricant density, kg/m 3 
dimensionless density, p/p Q 

density at atmospheric pressure, kg/m 3 
shear stress, N/m? 
dimensionless shear stress, t/E' 
shear stress ratio, t/t l 
limiting shear stress, N/m 2 
dimensionless limiting shear stress, t l /E' 

dimensionless Initial shear stress constant 


solid a 
solid b 

coordinate In direction of motion 
coordinate In direction of film thickness 
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THEORY 


The non-Newtonian approach will be to consider the flow conditions at the 
surfaces. The flow has two components, the flow due to velocity (Couette) and 
the flow due to the pressure gradient (Polseullle). 

In figure 2 we attempt to explain the velocity of the fluid for the five 
distinct zones that might exist In an elastohydrodynamlc conjunction when non- 
Newtonian effects of the lubricant are considered. In each zone the velocity 
of the top surface Is greater than that of the bottom surface. To better 
Indicate the difference between the zones, values of u ft and u b will be 
kept constant for each zone. In figure 2(a), zone 0, the normal Newtonian 
zone, the shear stresses at the surfaces are less than the limiting shear 
stress, and no slippage of the fluid at the surfaces occurs. In figure 2(a), 
zone 1, the Newtonian shear stress at the top surface Is larger than the 
11ml tiny shear stress and slippage occurs at the top (faster) surface. In 
figure 2(a), zone 3, the shear stresses at both the top and bottom surfaces 
are outside the limiting range (7 a < - and 7 b > 7^) and slippage occurs at 
both surfaces, but the slippage velocity Is less than the velocity at the sur- 
faces. In figure 2(b), zone 2, the Newtonian shear stress at the bottom sur- 
face Is larger than the limiting shear stress, and slippage occurs at the 
bottom (slower) surface. In figure 2(b), zone 4, the same sort of situation Is 
present as In zone 3 with the exception that the fluid slippage Is greater than 
the surface velocity. 

The relevant equations for the five distinct zones that can occur In an 
elastohydrodynamlc lubrication contact when non-Newtonian effects are con- 
sidered are developed next. Before we define these relevant equations for the 
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five zones, we need to define the dimensionless limiting sheer strength. For 
most materials the sheer strength verles linearly with pressure over wide 
pressure ranges, end there normally Is e certeln sheer strength even et zero 
pressure. It Is written es 

7j_ • Tq ♦ vP 0) 

where 

7q dimensionless Initial sheer strength 

y 11ml ting-shear-strength proportionality constant 

Jacobson (1970) found values of 7 q for e range of fluids to be 

between lxltT* 5 and 1x10“*. For the results given In that paper Tq 

was ass"«*ed tc be 9x10“^. Bair and Winer (1979) found y to be between 0.05 
and 0.10 for a complete range of natural and synthetic lubricating oils. 
Tevaarwerk (1976) used similar values of y In his traction utudles. In the 
present paper we consider a range of y between 0.04 and 0.10. 

Zone 0 

The Newtonian model gives 7 a < 7 L and 7 b < 7^. In this zone the 

rheological model Is a Newtonian fluid and no slippage occurs at either sur- 
face. The velocity distribution across the lubricant film for this situation 
Is depicted In figures 2(a) and (b), zone 0. The velocity, the velocity 
gradient across the film, and the shear stress at the surfaces can be written 
directly from Hamrock and Dowson (1981) as 
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equations (2) to (5) can be written In dimensionless form as 
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U-1-U*(1-2Z)+^V^Z(Z-1$ 


( 2 ) 

(3) 

(4) 

(5) 
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Since this Is zone 0, the Newtonian model gives 7 # < 7j_ ami 7 b < 7 L . 
The dimensionless mass flow per unit length and Its gradient set equal to 
zero give 

"■ 5dZ ' 7H ( 1 ‘iS;^) tu) 

dQ/dX - 0 gives 

&(?$)• " u (12) 

The second term on the right side of equation (11) Is the Polseullle or 
pressure term. For a Newtonian fluid equation (12) Is the familiar Reynolds 
equation when side leakage Is neglected. This equation will be referred to 
as the pressure equation since It will change for the respective zones. 


Zone 1 


The Newtonian model for zone 1 gives |7 a | > t l and |7 b | < 

The velocity distribution across the film that exists In this zone Is shown In 

figure 2(a), Slippage occurs at the top surface so that the velocity at the 
fluid-solid boundary Is denoted by u . This Implies that In this zone the 
velocity, the velocity gradient, and the shear stresses can be directly 

obtained from equations (2) to (5) with the substitution for u # 


with u,, the slip velocity. Making use of the non-Newtonian model discussed 

earlier leads to the observation that, since 7# > 7 l» then | Ta | 

Is set equal to t^. Solving for this slip velocity at surface a, we get 
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Making use of this equation, we can write the dimensionless velocity, 
velocity gradient, and shear stress at the bottom surface as 
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( 14 ) 


(15) 


■ ’L - ? Vjjls (16) 

The dimensionless mass flow rate per unit length and the pressure equation for 
this zone can be written as 


Q 





j£ 

dnU 



(17) 



The nonunity terms In parentheses In equation (17) are the Polseullle terms. 
Equation (18) Is referred to as the pressure equation but also can be viewed 
as the Reynolds equation for a non-Newtonian fluid when conditions In zone 1 
prevail. 


Zone 2 

The Newtonian model for zone 2 gives 7 a < 7 L and 7 b > 

The velocity distribution across the film that exists In this zone Is shown 
In figure 2(b). Slippage occurs at the bottom surface (slower moving 
surface) so that the velocity at the fluid-solid boundary Is denoted by 
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u fe . This Inpllss th«t In tin* zone tht velocity, the velocity gradient. 


and the shear stresses can be written directly f.aa equations (2) to (5) with 
the substitution for u^ with u^. Since the Newtonian Model gives 

*b > V the th€tr strtss at surface b Is set equal to t* l . Solving 

for the slip velocity at surface b, we get 



(19) 


Making use of this equation, we can write the dimensionless velocity, 
velocity gradient, and shear stress at the top surface as 



Note that, In order to be In zone 2, |T" a |< ^ rom •Nuatlon (22) we 
observe that this Is only possible If dP/dX Is negative. Likewise, In 
zone 1, |7j < 7 l . From equation (16) we observe that this Is only 
possible If dP/dX Is positive. 

The dimensionless mass flow rate per unit length and the appropriate 
pressure equation for this zone can be written as 
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7h ♦ t * 


IlH £ rirdp\ 

2nU Gbl™**/ 


(23) 
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■ 60 V? I 1 * 0*>{h<* H > - ■2JJ‘3 y( — = — ) 

The nonunity terms In the parentheses In equation (23) are the Polseullle terms. 


Zone 3 


The Newtonian model for zone 3 gives <r a > t L and equation (16) gives 
7 b < -7 l . Figure 2(a) shows the velocity distribution across the lubricant 
film. Slippage occurs at both the top and bottom surfaces* and the slip 
velocity is less than the velocity of the slow surface (surface b). In this 
zone the velocity* the velocity gradient* and the shear stresses can be written 
directly from equations (2) to (5) with the substitution for u. with 

“a and u b %. Since the Newtonian model gives t # > f L * the shear 

stress at the top surface, surface a, is set equal to 7 L and since 
equation (16) gives 7b - "H* the shear stress at surface b is set equal 
to -7 l . From these equations we find that u # « u b . The dimensionless 

velocity gradient can be written as 


H(Z* - Z)i L 


nU 


(25) 


* 

dU 

ST 


HT L (2Z - 1) 
nU 


(26) 


The dimensionless mass flow rate per unit length and the appropriate pressure 
equation can be written as 


Q 




(27) 


dP 

dX 



lL 

H 


(28) 


Th? nonunity terms In parentheses in equation (27) are the Polseullle terms. 
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Zone 4 
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The Newtonian model for zone 4 gives 7 b > 7^ and equation (22) gives 
7 a < -7^. Figure 2(b) shows the velocity distribution across the lubricant 

film. Slippage occurs at both the top and bottom surfaces. The difference 
between zone 3 and zone 4 Is that In zone 3 the slip velocity Is less than the 
slower surface velocity and In zone 4 the slip velocity Is greater than the 
faster surface velocity. In zone 4 the velocity, the velocity gradient, and 
shear stresses can be written directly from equations (2) to (5) with 

the substitution for u with u. and u K with u K . Since equation (22) 

a a D d 

gives £ -7^, the shear stress at surface a Is set equal to -7^; and 
since the Newtonian model gives 7 b > 7 L , the shear stress at surface b Is set 
equal to 7 L . The dimensionless velocity, velocity gradient, mass flow, and 
pressure can be directly written from the equations developed In zone 3 
(eqs. (25) to (28)) If ~ L Is substituted for 7 L . Note that In 
order to be In zone 3 from eauatlon (28), the pressure gradient must be 
positive; and In order to be in zone 4, the pressure gradient must be 
negative. 


Oensity, Viscosity, and Film Shape Equations 


The eauations that define velocity, velocity gradient, shear stress, flow, 
and pressure for the five zones have been defined. Before proceeding, however, 
the dimensionless density, viscosity, and film shape need to be expressed. For 
a comparable change In pressure the density change Is small as compared with 
the viscosity change. However, very high pressure exists In elastohydrodynamlc 
films, and the liquid can no longer be considered as an Incompressible medium. 
From Dowson and Higglnson (1966) the dimensionless density can be written as 



1 + 


0.6E'P 

I ♦T.TET 


(29) 
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where E 1 Is expressed In glgepascals. 

The effects of pressure on viscosity can be written as 

. i. . «*P . « 6P ORIGINAL PAGE 13 (30) 

OF POOR QUALITY 

where 

a pressure-viscosity coefficient of lubricant, m^/N 
rig viscosity at atmospheric conditions, N s/m 2 

The film shape equation can be written from Hamrock and Jacobson (1982) 

In dimensionless form as 

H k -** H o + i[r(r) + 4^ p i D j] (31) 

where 

j - k - 1 + 1 

D » b[fX - B)1n(l - B) 2 - (1 + B)ln(l ♦ B) 2 ♦ 4B(1 - In b)] 
and 

b semi width of the Hertzian contact, R V8W/w 
B 1/n 

n number of nodes In semi axis of contact 

The last term In equation (31) represents the elastic deformation at a 
point x due to the contribution of various rectangular areas of uniform 
pressure in the conjunction. 

Having defined density, viscosity, and film shape, we can return to the 
solution of the mass flow and pressure equations for the five zones. The first 
step Is to rewrite these equations by finite difference approximations which 
rely on the fact that a function can be represented with sufficient accuracy 
over a small range by a quadratic expression. Standard finite central- 
difference representation was used, and a procedure similar to that given In 
Hamrock and Jacobson (1982) was followed. Figure 3 shows the flow chart of the 
pressure loop In the computer program used to evaluate the mass flow and pres- 
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sure from the Inlet to the outlet. In this figure the procedure used In 
evaluating zone 0. the Newtonian model, was quite similar to that used In 
Hamrock and Jacobson (1982). A similar procedure was used In evaluating the 
other zones. 

Having defined the pressure from the Inlet to the outlet for the non- 
Newtonian fluid model, we can evaluate the force components, the shear forces, 
and the coefficient of friction. The relevant expressions are presented In 
Hamrock and Jacobson (1982) and are not repeated here. 


Cr 
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RESULTS 

Olnonslonless Grouping 


From the variables of the numerical analysis the following dimensionless 
groups can be defined: 


(1) Dimensionless film thickness: H - 

w z 

(2) Dimensionless load parameter: W ■ -£t r 

n 0 u s 

(3) Dimensionless speed parameter: U ■ 

(4) Dimensionless materials parameter: G « oE* 


(5) Dimensionless sliding velocity: U* « ± ,, b 

u s u a u b 

(6) Limiting-shear-strength proportionality constant: y 

The first four groups were used for the Newtonian rluid analysis In Hamrock and 
Jacobson (1982). the dimensionless film thickness for a rectangular contact with 
a non-Newtonian fluid can be written as 


H min * h m in, N f ( u *-’' u - w - G > 


(32) 
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The dimensionless minimum film thickness for • Newton 1 §n fluid as obtained from 
ilamrock and Jacobson (1982) can be written as 

H m1n,N " 3,07 U°* 71 G°* 5 V 0,11 (33) 

In the present analysis the dimensionless parameters U* and y were varied 
and the effect on the minimum film thickness was studied. 

Pressure and Film Profiles 

Representative variations of dimensionless pressure, film thickness, shear 
stress distribution, and mass flow are shown In figures 4 to 6. In these 
figures the Inlet region Is to the left and the outlet region Is to the right. 
Figure 4 shows the variation of dimensionless pressure and film thickness on 
the X axis for four dimensionless sliding velocities. The Hertzian pressure Is 
also shown In this figure. The characteristic pressure spike Is clearly 
evident for each dimensionless sliding velocity, but the spike diminishes as 
the sliding velocity increases. Also clearly evident for each sliding velocity 
is the parallel film shape through the central part of the contact, with a 
minimum film thickness occurring near the outlet of the contact. As the dimen- 
sionless sliding velocity Increases, the central parallel film thickness 
decreases more than the minimum film thickness decreases. As a result the nip 
decreases with Increasing sliding velocity, where the nip Is that portion of 
the film shape from the tip of the pressure spike to the outlet of the 
conjunction. 


Variation of Shear Stress 

lhe shear stress ratio Is defined as the shear stress divided by the 
limiting shear stress. Variations of the shear stress ratio at the top 
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{t a ) and bottom (t b ) surfaces on the X axis are shown In figure S for the 

four dimensionless sliding velocities given In figure 4. Recall that the top 
surface (surface a) Is moving with a greater velocity than the bottom surface 
(surface b). For the condition of no sliding (U* ■ 0) zone 0, the Newtonian 
model. Is valid for the entire conjunction, and the shear stress ratio of the 
top surface t . Is symmetrical to the shear stress ratio of the bottom 

fit 

surface about the zero axis. From equations (9) and (10) this effect 
Is clearly evident. In these equations the first term on the right Is set 
equal to zero since U* - 0, and the only remaining effect Is the pressure 
gradient. There Is a difference In sign for the different surfaces. For 
sliding velocities U* of 0.01 and 0.02 zones 1 and 2 appear at the tip of the 
spike and zone 0, the Newtonian model, occurs elsewhere. For U* ■ 0.04 the 
zone alloration going from the inlet to the outlet (left to right) Is zone 0 to 
zone 1, to zone 2, back to zone 0; then In the spike zones 1 and 2 reoccur. 
Zones 3 and 4 do not occur in the final results for figure 5, but these zones 
are used to derive a converged solution. For the higher load cases, however, 
zones 3 and 4 do appear in the final converged solution. 

The shear stress ratio in the Inlet up to the beginning of the Hertzian 
contact Is quite similar for all four sliding velocities. From the beginning 
of the Hertzian contact to the tip of the pressure spike there Is substantial 
change In the shear stress ratio from U* » 0 to U* • 0.04. As the dimen- 
sionless sliding velocity increases, the shear stress ratio Increases until 
for U* » 0.04 there Is a large region where the shear stress Is equal to the 
limiting shear stress. From figure 4 we found that the pressure spike de- 
creases with increasing sliding velocity. The result of this Is clearly 
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evident In figure 5, where fro* the tip of the pressure spike to the exit the 
sheer stress ratio at the top surface t a becomes less negative with In- 
creasing dimensionless sliding velocity. The effect of using the non-Newtonian 
model Is quite significant for the largest dimensionless sliding velocity In 
figure 5. 


Flow Results 

A typical variation of dimensionless flow and the Polseullle flow terms 
on the X axis Is shown In figure 6. Since there was very little difference 
for the various sliding velocities, a typical result Is shown. Equations (11), 
(17), (23), and (27) define the dimensionless mass flow rate for the five zones 
considered. The Polseullle flow term In these equations Is the term or terms 
other than the unity In the parentheses on the right side of these equations. 
From figure 6 we find the flow to be constant throughout the contact. Great 
care was taken to assure that this was maintained for all the results pre- 
sented. As In the authors' earlier publication (Hamrock and Jacobson, 1982) 
slight adjustments In pressure profile were necessary In the Inlet region to 
assure a constant flow. The Polseullle term approaches unity at the Inlet, Is 
near zero from the beginning of the Hertzian contact to the location of the 
pressure spike, and Is negative from the pressure spike to the outlet. 

Effect on Coefficient of Friction 

Figure 7 shows the effect of dimensionless sliding velocity U* on the 
coefficient of friction » for three dimensionless loads. For higher loads 
the asymptotic value of the coefficient of friction Is reached at lower dimen- 
sionless sliding velocities. The dimensionless sliding velocity given in this 
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figure can be equated to the conventional sliding velocity In percent by 
multiplying the value of U* by 200. For example U* ■ 0.03 corresponds to 
a conventional 6 percent sliding speed. The results for a Newtonian fluid 
are shown In this figure with dashed lines. The asymptotic value of the co- 
efficient of friction Is between 0.07 and 0.08 for a 11 ml ting-shear-strength 
proportionality constant y of 0.07. 

Influence of Sliding Velocity 

Cases 1 to 6 of table I show the effect of six values of dimensionless 
sliding velocity U* on dimensionless minimum film thickness for constant 
values of the other dimensionless parameters. From these results It Is ob- 
served that as the sliding velocity Is Increased the minimum film thickness 
decreases. However, the Influence of sliding velocity was found to be not 
simply a function of minimum film thickness but had to be coupled with the 
dimensionless speed, load, and materials parameters. Table I also shows these 
results. The resulting relationship Is 

^ min . - « 1 - 4.07 x10" 9 (U*) 0,60 U 0,23 (WG 2 ) 3 ’ 85 (34) 

^min,N 

The Influence of the maximum pressure Is represented In this equation by the 
grouping WG^. The coefficient of determination e was also calculated for 
these results. The value of c reflects the fit of the data to the resulting 
equation: unity representing a perfect fit, and zero the worse possible fit. 
The coefficient of determination for these results was 0.9990, which Is 
excellent. 
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Influence of Limiting Shear Strength 

Cases 9 and 10 of table I show the Influence of the llmltlng-shear- 
strength proportionality constant y on dimensionless film thickness for 
constant values of the other dimensionless parameters. From these three cases 
covering the complete range of natural and synthetic lubricating fluids, the 
following relationship was found: 

^12- « 1 ♦ 2.06{y - 0.07) (35) 

^m1n,N 

The coefficient of determination was found to be 0.9995, which Is excellent. 

Minimum Film Thickness 

The proportionality equations (34) and (35) have established how the non- 
Newtonian fluid rheological model affects the minimum film thickness for low 
sliding speeds. At high sliding speeds the oil film thickness Is mainly 
governed by the velocity of the slower surface. The reason for this Is that 
the shear strength Is reached at the faster surface and the oil velocity Is 
given by the slower surface. This phenomenon occurs In the Inlet of the con- 
junction. The hlgh-sli ding-speed cases have not been numerically computed 
because of numerical Instability. For the film thickness equation to accommo- 
date these hlgh-sli ding-speed situations, the film thickness equation had to be 
modified as Indicated below. 

H m1n - N (exp [-4.09 x 10 _9 (U*) O - 6O U°* 23 (WG 2 ) 3 ' 85 ♦ 2.06(r - 0.07)] 

♦ u4 ),71 (l - U*) 0 - 71 (36) 
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This dimensionless mlnlmum-1 sothermal-f Urn thickness formula Is for .ully 
flooded rectangular lubricated contacts Incorporating non-Newtonian rheological 
effects. Asymptotically this equation gives the minimum film thickness at high 
pressures (>1 GPa as equal to zero If one of the surfaces Is stationary. 

Table I gives the fourteen cases used In evaluating equation (36). In 
this table H^, the minimum film thickness obtained from equation (36) 9 
corresponds to the minimum film thickness obtained from the non-Newtonian 
elastohydrodynamlc lubrication theory developed earlier. The mlnlmum- 
fllm-thlckness equation obtained from Hamrock and Jacobson (1982) for the 
Newtonian elastohydrodynamlcally lubricated rectangular contact Is denoted 
by and Is given In equation (33). The ratios H vn1n /T^ n1nv N 

and flmin/^nln.N are a1so 9 iv * n 1° table !• The percentage difference 
between and Is expressed by Vj» which Is defined as 

Vj . ( ) loo (37) 

In table I the values of are within *1.0 percent. 

ORIGINAL PAGE IS 
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The values of the coefficient of friction for the 14 cases studied are 
given In table I. Making use of these results, we can write an approximate 
formula for the coefficient of friction as 

3.32 

. 0.67x 10" 6 (U*)°* 81 U°’ 26 (WG 2 ) ’ when ^ < 0.8 * (38) 

If the coefflcent of friction ^ Is greater than 0.8 ft the following 
approximate expression should be used: 

w . 0.80 v * 0.021 tanh - 0.8o] when ^ > 0.8 y (39) 
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This coefficient of friction Is mainly determined by the sheer strength of 

the lubricant. The percentage difference between u and u Is expressed 
by V*. where 

V 



In table II the values of Vj are within *9 percent. Equation (39) was 
derived from eight calculations of the coefficient of friction with sliding 
velocities up to 20 percent and different y values. In these calculations 
the pressure distributions from the lower sliding velocities were used. 

CONCLUSIONS 

A procedure for the numerical solution of the complete elastohydrodynamlc 
lubrication of rectangular contacts Incorporating a non-Newtonian rheological 
model Is outlined. The approach uses a Newtonian model as long as the shear 
stress Is less than the limiting shear stress. If the shear stress exceeds the 
limiting shear stress, the shear stress Is set equal to the limiting shear 
stress. The limiting shear stress Is expressed as a semlemplrlcal linear 
function of pressure. The numerical solution therefore requires the coupled 
solution of the pressure, film shape, and fluid rheology equations from the 
Inlet to the outlet without making any assumptions other than neglecting side 
leakage. 

By using the procedures outline In the analysis the Influence of the 
dimensionless speed U, load W, and materials G parameters, dimensionless 
sliding velocity U+, and limiting shear-strength proportionality constant y 
on minimum film thickness H m1n has been Investigated. Fourteen cases were 
used to generate the minimum film-thickness equation for an elastohydro- 
dynamlcally lubricated rectangular contact Incorporating a non-Newtonian 
rheological model: 
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3.85 

*»1n “ ^1n t N {•xpH.07xlO- 9 (U*) 0 - 60 U°- 23 («6 2 ) * ♦ 2.06(y - 0.07)3 

♦ (1 - 

The mlnlmum-fllm-thlckness equation obtained for an elastohydrodynamlcally 
lubricated rectangular contact Incorporating a Newtonian fluid rheology model 
'was developed by the authors In an earlier publication (Hamrock and Jacobson, 
1982) where 


mln.N 


3.07 u0.71g0.57 w -0.11 


Besides the dimensionless film thickness formula, formulas for the coefficient 
of friction were developed. 

U] • 0.67k10 _6 (U»)°' 81 U O * 26 ( HG 2 ) for ITj < 0.8 t 

r-i 1 

v ■ 0.080 y ♦ 0.021 tanh I—— - 0.80 I for > 0.8 y 

Computer plots are presented that Indicate In detail pressure distribution, 
film shape, shear stresses at the surface, and mass flow throughout the con- 
junction. The characteristic pressure spike Is clearly evident for each of the 
sliding velocities, but the spike diminishes as the sliding velocity Increases. 
Also clearly evident In the computer plots Is the parallel film shape through 
the central part of the contact, with a minimum occurring near the outlet of 
the contact. The central parallel film thickness decreases more than the 
minimum film thickness decreases with Increasing sliding velocity. The result 
of this Is that the nip decreases with Increasing sliding velocity, where the 
nip Is that portion of the film shape from the location of the pressure spike 
to the outlet of the conjunction. The computer plots of the shear stress at 
the surfaces Indicate that as the sliding velocity Increases these stresses 
approach or equal the limiting shear stress. A simple computer plot of the 
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flow shows that It Is constant throughout the conjunction. The effect of 
sliding velocity on coefficient of friction for various loads Indicates that 
as the load Increases, the limiting value of the coefficient of frlcton Is 
reached for a much lower value of sliding velocity. 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135, December 23, 1982 
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TMU I. - EFFECT OF DINERS IOMLESS MJMMETEJIS ON 


Cm* 


Oiaenslcn- 
1m loa* 
par anetcr, 
M 


Oh 
l«« tpw4 
par aaeter, 

U 


Olnentlonless 

naterlals 

HTMtir, 


Dtamfo*- 
l«I sliding 
velocity, 
O* 


Uniting- 

shear- 

strengtli 


tlonalltjr 

constant, 

r 


Oil 

less nlnlMi 
file tMctress, 
Nrin 


■Ssln 


Knt*,R 


1 12.0479*10-* 


1.0000*10-11 


sooo 


12 

13 

u 


]4.oo*4 

13.0000 

fo.oooo 


5.5975 
5.5975 
t.a 
1. 


To 


0.07 


.020 

.030 

.040 


34*1.1 

3491.1 

SOOO 

sooo 


.040 

.020 

.020 


.014 


.04 

.10 

.07 


.004 


19.524*10-* 

19.300 

19.294 

19.240 

19.974 

19.497 

20.327 

20.133 
19.449 

20.134 
52.502 
51.049 
19.044 
19.442 


1.000 


.972 

.947 

1.000 

.990 

.944 

1.031 

1.000 

.973 

1.000 
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TABLE II. - EFFECT OF OIMENSIONLESS PARAMETERS ON 
DIMENSIONLESS CENTER OF PRESSURE 


Case 

Dimension- 
less center 
of pressure, 
X C p 

Dimension- 
less center 
of pressure 
from least- 
squares fit, 

x cp 

Difference, 

V 2* * 
percent 

Coefficient 
of friction, 
y 

Coefficient 
of friction 
from least- 
squares fit, 
y 

Difference, 

V3 ’ ♦ 
percent 

1 

2 

-0.1487 
- .1480 

-0.1487 
- .1469 

0 

- .71 

0.00054 

.01227 

0.01250 

1.87 

3 

- .1430 

- .1452 

1.52 

.02015 

.02191 

8.73 

4 

- .1384 

- .k<’ 

2.35 

.03613 

.03842 

6.34 

5 

- .1361 

- .l'oul 

1.49 

.05309 

.05336 

.51 

6 

- .1346 

- .1346 

0 

.06316 

.06403 

1.38 

7 

8 

- .1912 

- .1845 

- .1912 

- .1730 

0 

- 9.48 

.00081 

.03479 

.03210 

-7.73 

9 

- .1360 

- .1417 

4.16 

.03618 

.03534 

-2.32 

10 

- .1383 

- .1417 

2.46 

.03592 

.03842 

6.96 

11 

12 

- .1810 
- .1841 

- .1810 
- .1753 

0 

-4.79 

.00202 

.04678 

.04655 

- .49 

13 

14 

- .0971 

- .0973 

- .0971 

- .0960 

0 

-1.34 

.00087 

.04744 

.04441 

-6.39 
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Figure 3. - Flow diagram. 
















(i) u * ■ o. o>> u* - aoi. 

lc)U**(X02. (dHJ* ■ OL04. 

Figure 4 - Prtssur* distributions and film shapas tor olffarant sliding velocities U*. Dimansionlass load panmatar W, 2,048x10**; 
dimensionless spaad paramatar U, UT 11 ; dimensionless materials parameter G, 5000; llmltlng-stieer-strangth proportionality 
constant y, 0.07. 
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NODES ON X-AXIS 

<•» u* • o. «>) u* - aoi. 

(c)u*-o,(». Miu**aoi 

Figure S. - Shaar stress distributions and lubricating waiter different sliding vttocltias U\ Dlmanslonltss load parematar 
w, 2.0U8xl(T 5 ; dlmanslonlass spaad parematar U, 10' 11 dlmanslonlass matt rials parematar C, MOO; llmltlng-shaar- 
strength proportionality constant y, a OF. 
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DIMENSIONLESS 
LOAD PARAMETER, 

/ 3. OOOOKlO' 5 / 


f 2. 0478x10 


I63B2X10" 5 . 


NEWTONIAN OIL 


0 .01 .02 .03 . 04 . 0$ .06 .07 .OB .09 .10 

DIMENSIONLESS SLIDING VELOCITY, U* 


0 24 68 1012 14 161820 

SLIP, PERCENT 

Figure 7. - Coefficient of friction as a function of load and dimen- 
sionless sliding velocity. Dimensionless speed parameter U, 10 , 

dimensionless materials parameter G, 5000; limiting-shear* 
strength proportionality constant y.0.07. 


